Dephasing and the steady state in quantum many-particle systems 
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We discuss relaxation in bosonic and fermionic many-particle systems. For integrable systems, the 
time evolution can cause a dephasing effect, leading for finite subsystems to certain steady states. 
We give an explicit derivation of those steady subsystem states and devise sufficient prerequisites 
for the dephasing to take place. We also find simple scenarios, in which dephasing is ineffective and 
discuss the dependence on dimensionality and criticality. It follows further that, after a quench of 
system parameters, bipartite entanglement entropy will become extensive. This provides a way of 
creating strong entanglement in a controlled fashion. 

PACS numbers: 05.70.Ln, 05.30.Ch, 02.30.Ik, 03.67.Bg 



In equilibrium statistical physics we usually work with 
canonical ensembles, characterized by density matrices 
Q cx 6^^''"'='^'' that can be derived by maximizing the 
entropy under the constraint of fixed expectation values 
(Ofc) for a few (macroscopic) observables like energy or 
particle number However, it is in general unclear 
whether a given system in a certain initial state will 
evolve to any steady state at all and if so, whether that 
state is indeed one of the canonical ensembles. Recent ex- 
periments, especially with ultracold gases, have revived 
the interest in this important topic of nonequilibrium 
physics. In such experiments with integrable systems, 
absence of thermalization was observed (e.g. in z|). We 
will deduce analytically under what circumstances inte- 
grable systems relax to non-canonical steady states. 

It was conjectured in Q] that the time evolution in an 
integrable system with conserved observables Ik should 
lead to the corresponding maximum entropy ensemble [IJ 

ed = |e^^-"'=^"^ (1) 

where the ak are determined by the initial state. So 
far, the conjecture was discussed by analyzing specific 
observables for specific systems 0, i, S i]. In 01 it 
was found that certain observables do not relax to the 
ones predicted by the ensemble gd as conjectured in 
So here the focus will be not on observables but states 
themselves. 

We point out that it is essential, for the conjecture to 
hold, to restrict oneself to measurements in a finite sub- 
system, i.e. dl]) should not be interpreted as the steady 
state of the full system (see a first discussion in Q). A 
general proof is given, showing explicitly how state relax- 
ation in such subsystems can occur due to dephasing. We 
devise the necessary prerequisites, discuss the relaxation 
speed, and give simple examples and counter examples 
for dephasing. Our results allow for a simple interpreta- 
tion of the discrepancies between Q and 0] ■ On general 
grounds it follows further that through dephasing the bi- 
partite entanglement entropy in a pure state will become 
extensive. This may be of interest for quantum compu- 
tation applications where entanglement is needed as a 



resource. We comment on implications for the ability to 
simulate such systems on classical computers. We will 
first present the case of free systems, then point out how 
this generalizes to Bethe Ansatz integrable models and 
close with a short discussion of nonintegrable systems. 

Dephasing for quadratic Hamiltonians.- First, we ad- 
dress the situation of quadratic Hamiltonians (for t > 0), 

H = E^J [alV^Ja^ + i(a|VKya] + h.c.)] + const. , (2) 

where are bosonic or fermionic (C :— ±1) ladder oper- 
ators, [Cj , aj]-c = Sij, and = (oi, 02, . . . ). This covers 
also lattice regularized free field theories. The Hamilto- 
nian is diagonalized by a linear canonical transformation 

O-c/^t), H^iv^ysv ^Ek^kviVk^ (3) 

where e is the diagonal matrix of one-particle ener- 
gies. For quadratic systems, the groundstate, thermal or 
evolved states are all Gaussian. Assuming such a Gaus- 
sian initial state gt=a, due to the Wick theorem, the state 
is always fully characterized by the one-particle Green's 
functions Q (superscripts a, rj indicate the chosen basis) 

g^giG), g^^{{:,)-{'^^)'),^Ug^W. (4) 

The initial state go might e.g. be the groundstate of a 
different quadratic Hamiltonian ( "quench" of system pa- 
rameters at t = 0). In the following, we will consider 
bipartitions of the full system into a subsystem O and its 
environment fi"*", call Oq the projection onto and the 
volume V :— Volf2 -I- Voin-*-. In the situation where O-*- 
is approaching the thermodynamic limit, the quantum 
number labels k become a set of d continuous labels and 
one discrete label k (fc, s) (e.g. a d-dimensional mo- 
mentum vector and a band index, Ek = £fcs denoting the 
dispersion relation) with k gT, VolF finite, and the den- 
sity of states ps oc 1/V. Again, due to the Wick theorem, 
reduced density matrices Tr^i g of quadratic systems are 
functions of the one-particle subsystem density matrices 
G, TTn± g = g{G), where G is defined by = QnG^^n- 
Theorem. Let the Green's function Qd be defined by 

[Qd]kk' := 4fe'KLo]fcfc ' (5) 
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i.e., in the eigenbasis representation, Qd is the diagonal 
part of Gt=o- Under preconditions (a-c), stated after the 
proof, the steady t — > oo state of O is then given by 

lim Trf2_L Qt = TTn± Qd with Tr^-L Qd = g{Gd) , (6) 

t — >oo 

i.e. subsystem states relax to the reduced density matri- 
ces of Qd = Qd, the maximum entropy ensemble ([l} with 

h^^W and a.. = ln(C/(l-[t;2]^,i)). (7) 

Proof. To compare the two reduced density matrices 
in proposition ([6]), we compare at first the corresponding 
subsystem Green's functions hmj^oo Gf , Gd- They are 



= Qn{Ug'lU^)Qn , and 



(8) 



where = e^*/*'*, = e^*/^'^, £ ;= _^), and thus 

G? = ^Efcfc'e(^'=-^'=')*M/fcfc', with (9) 
/fefc' V • enJ7fc [^o^lfefe' Ul,Qn, [Uk]^ C/.fc . (10) 

Here, fkk' is a matrix valued function of k and k' , whose 
matrix indices label points in Vl. Comparing G° to G^, we 
see that ^ is true, if the nondiagonal contributions k ^ 
k' to G^ vanish for t ^ cxd. The summation over /c, for a 
fixed Afc = k — k' ^ 0, corresponds for increasing i to a 
Fourier transform with respect to ever higher frequencies 
which may vanish due to phase averaging (hence we call 
the effect "dephasing"): To see this, let us rewrite G°, 
(HD, with (fc, k') (fc, s,k' = k + Afc, s') as [13 



G? 



E — 



Afc 



d''fce-''^^''=('^)*/A'fc(fc) 



(11) 



with phase function ip and group velocity difference g, 
<^Afe(fc):=^^^^%^, 9ll(fe):=afeV'll(fe). (12) 



In the following, we omit indices s, s', and Afc and con- 
sider always nondiagonal contributions to pT|) . i.e. s ^ s' 
or Afc 7^ 0. If € C^, g finite on supp(/) C F, and if the 
matrix elements of / are integrable [16|, the integral 
of (fTT|) vanishes for t — > oo (Riemann-Lebesgue Lemma) . 
This is due to the fast oscillation of the Fourier kernel 
g-ig (fc-fco)* the vicinity of every point fcg. Further, if 
/ is bounded in a vicinity B of fco, the contribution of 
B to the integral is of 0{yo\B /t). But what happens at 
points fco where g vanishes? If such a zero of g, is isolated 
(i.e. if the Hesse matrix of ip at fcg is invertible) and if / is 
smooth in a vicinity of fco, its contribution to (|lip is still 
vanishing (follows from the Morse lemma and a Gaussian 
integral) . We can also treat the case where the Hesse ma- 
trix := mimVvWo W^hW has only n (1 < n < 
d) nonzero eigenvalues {/ii}i<i<n- With the transforma- 
tions Mi-.k^ {q', q") := q = W ■ {k - ko) {q' e M", 



q" g Rd-n) and 7W2 : g -> (Q', Q") = {h'^^'^q', q"), the 
contribution from a vicinity B of fco becomes (see Fig. [1]) 

/gd'^fce-^^W*/(fc) d'^ge-^*2:Q?+0(Q^)/2(Q) 

- /dPe-"-PF(P), (13) 

where B2 = M2{Mi{B)), h ^ f ° o A^-^, and 

^(^) ^= iTP /B„|Q'|=ypd'^"'Q/2(Q)- Ii^tegral m has 
again the form of a Fourier transform and vanishes for 
t ^ 00, a F G £^(^2). If F is bounded in the vicinity of 
P = 0, the integral vanishes as 0{Vo\B/t). Those con- 
ditions on F have to be interpreted as stricter conditions 
on / when approaching zeros of g, appropriate to still 
guarantee the dephasing. Zeros of g with multiplicity 
i — 1 > {Sj^ = 0) can be treated in a general fashion 
only for Id by substituting Q = q^. In higher dimensions 
they are a lot more complicated; see e.g. [§] for 2d. 

In the cases discussed so far, we have assumed that tp 
is continuously differentiable {(p S C^). For the sake of 
brevity we will consider now only the scenario y'(fe) = 
ipo + \kf + 0(/fc^+i). This is for ^ = 1 nondifferentiable. 
The cases £ = 1 and £ = 2 cover very typical examples 
(for example magnons, phonons or free particles). The 
contribution of a (small) sphere of radius K around fco = 
to the integral in (fTT|) is with the substitution P = k^ 



dPe'^f°+P^*F{P), F{P) 



P- 



d'^-;fc/(fc). 

|fc|=P7 

(14) 

00 if F e Cj^^ and if 



This contribution vanishes for t 
F is bounded, the contribution is of 0{K'^/t). 

For a few general, physically relevant scenarios, we 
have established sufficient conditions under which (every 
matrix element of) the integral in (fTT|) goes to zero as 
t — > 00. As Vol F is finite and as we have a finite number 
of bands s, those conditions guarantee hence that all non- 
diagonal contributions to Gf vanish. Due to Wick's the- 
orem, expectation values of arbitrary observables O on f2 
are given by polynomials in [G°]y . Nondiagonal contri- 
butions to (O) will consequently also vanish if we restrict 
to finite subsystem sizes. From the convergence G° 
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FIG. 1: Sketch of coordinate transformation and integration 
path in (|13[) for a non-isolated zero of g with d = 2, n = 1. 
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follows then the proposition ([6]) . Finally, ^ follows from 

= l/(e"'' - C) and avbk)s. = K'lfcfc " 1- □ 

Preconditions. During the proof, we collected the fol- 
lowing prerequisites, sufficient to guarantee convergence 
to the steady state dH]): (a) Vol 17 is finite and V ^ oo. 
(b) The parameterization fc ^ (fe G F, s) of the quantum 
numbers is possible with a finite Vol F and a finite num- 
ber of bands s. (c.l) Lp g C^, g finite on supp/ C F, and 
f £ . If this is not given for the vicinity of a point fcoj 
we require for such a point (c.2) if g{ko) = 0, then the 
Hesse matrix 9j(p\ka exists and is nonzero, and F E L^, 
ov (c.3) if, at fco = 0, ip(k) = ipo + \k\^ + ©(fc^+i) then 
F g Of course, (c.3) can be generalized to the case 
\k\^ — > |A • (fc — ko)Y with some nonzero matrix A. 

With those conditions, all nondiagonal contributions 
to the subsystem Green's function matrix G° vanish for 
t ^ oo and dephasing to the steady state ensemble is 
effective, ([6|). If /, i^, and F are, in the correspond- 
ing situations, bounded instead of only integrable, 
nondiagonal contributions to decay (more quickly) as 
0{l/t). Below, we give illustrative examples. Among 
those are simple scenarios where some of the prerequi- 
sites are violated and dephasing does in fact not occur. 

Examples and counter-examples for dephasing.- At the 
end of the proof, a reason for requiring (a) was given. 
A simple counter-example consists in violating (a) with 
fi-L = and measuring C (77^, 77^ = e"*'^"'*[C/^]fcfc'-4fc', 
i.e. measurements in infinite subsystems, can reveal the 
phases and nondiagonal contributions ("rephasing"). 

If g has zeros or if / has divergences, dephasing proper- 
ties are dominated by the vicinities of such points. Thus, 
we illustrate (c)hy considering the paradigmatic scenario 
tf{k) ~ (/?o + |fc|^, f{k) ^ l/fc™ near fc = (for some fixed 
s, s', Afc and i.j E ft). The integral in pT|) is then 



,iqt 



_ m+e-d 



(15) 



Hence this (nondiagonal) contribution to G^, for t — > 00, 
does not vanish if x > 1, vanishes as l/t^~^ if < x < 
and (at least) as 1/f if x < 0; see Fig.O In this scenario, 
both (c.2) and (c.3) apply with F,Fq(: P^^ and are not 
only sufficient but also necessary. 

Our first explicit example is {ai = Cx, h = 1) 



x+l 



(16) 



the dimerized fermionic tight-binding model, where 
modes and Ck+TT are coupled and the dispersion rela- 
tion is e/c± — ±2\/cos^ fc + 7^ sin^ fc, i.e. gapless if 7 = 0. 
The eigenmodes are labeled rik±. We evolve the ground- 
state for a certain dimerization 70 with a different value 
7 7^ 7o. Skipping details of the calculation, we note that, 
using (fTTj) . the nondiagonal contributions (c^c^/)"^ from 
iVk+vi-) and {r]f^_ril^) to {Cxc\,)t can be written as 



t^"'^-/^dfc/(fc).^°^('^(^)*)' °ddx-.'^ ^^^^ 



> a) ^ = 1 1 

^ ^^^---f oc 


dephasing 
if m < 1, 
dephasing 0(1/4) 
[ if m < 0. 
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FIG. 2: Pliase function tp, sin(<^(fc)t), and the nondiagonal 
contribution / to G" in the paradigmatic case (|15p . d= \. 
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FIG. 3: a) Dimerized fermionic Id tight-binding model (|16p . 
Dispersion relation ek± = ^^ip{k)/2 for 7 = 1, |, i,0 and the 
nondiagonal contribution f{k) to {cicj)t according to Eq. (|17|) 
for quenches 7=i— >0, |^|,0— >| (top to bottom), b) 
The d- dimensional harmonic lattice model. Dispersion rela- 
tion Ek = — f{k)/2 for Lo = |, |,0 with 7 = 1, and the non- 
diagonal contribution f(k) to {hxhl^,)t according to Eq. (|19p 
for quenches tJ = |— >0,0— »|,|-^| (top to bottom). 



Fig.[3K displays ek± = ^ip{k)/2 and / for [x, x') ~ (1,0) 
and three different quenches with 7 7^ 1. The zeros of 
g are k = 0,±-| and have each, in the notation of the 
paradigmatic situation (|15p . one of the characteristics 
(l,m) = (2, -2), (1,0), (2,-1). Hence, x < and de- 
phasing of 0{\/t) is guaranteed. The same is given for 
all other {x,x') and has also been checked numerically. 
However, if we switch to 7 = 1, is const. = — 4 Vfc and 
no dephasing can occur ~ we have uncoupled dimers. 

As a second explicit example we choose the harmonic 
lattice model in d dimensions (a^ = 5^ , fi. = 1). Contrary 
to the first example, it will not dcphase in all cases. 

H = \ Y.r\Pl + ^'0' + 7 Etl {Qr - Q.+eJ^] • (18) 



i sin((y9(fc)i), even x — x' 



The dispersion relation is £fc = y -I- 47 Ylii=\ ^'^'^ ki/2, 
i.e. gapless for lo = 0. With the bosonic operators br ■= 
{Qr + iPr)/V^, H is brought to the form ^ and is hence 
amenable to the theorem. Using (lll|) we arrive at 

{Kbl)f = 11^ d'^kfik) cos(2efet), / = e^Mr-r')j^ 

where we switch at i = the oscillator frequency u;' — s- w 
and is the dispersion relation before that quench. The 
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dephasing properties are dominated by the vicinity of 
fe = 0. If one switches between two noncritical values 
Lo' ,uj > 0, its characteristic in terms of p3|) is (^,m) — 
(2,0), i.e. X = and hence dephasing of 0{l/^/t) for 
d= \ and 0(l/t) for d > 1. For = 0, w > 0, one has 
(Z,to) = (2, 1) and consequently no dephasing for d — 1, 
and dephasing of 0(1/ v^) [0{l/t)] for d = 2 [d = i]. 
For > 0, a; = 0, one has (/,to) = (1,2) and hence no 
dephasing for d = 1,2, and dephasing oiO{l/t) iov d — i; 
Fig. [3b. This was confirmed numerically for d = 1, 2. 

As a last example, consider free hard-core bosons in a 
Id box. The system is prepared in the groundstate for 
a box of size L which is switched to L > L at t = 0, 
3 . The Jordan- Wigner transformation yields a model of 
free fermions. The transformation between one-particle 
eigenstates before and after the quench {\q) and \k)) is 



Vkq = {k\q) 



^i(k-q)/2 sin(£(fc-g)/2) 
sin((fe-g)/2) 



'LL 



Vk-q. (20) 



The weight of Vaa: is concentrated in the interval |AA:| < 
2tt / L. With the Fermi momentum qp, the initial Green's 
fimction [^g]gq' = 6qqi {\—6{qp — \q\)) is diagonal in the \q) 
basis and [t/o]fcfc' = \VQ'^V^kk', which appears in ([10]), is 
hence also concentrated in |fc— k'\ < 1/ L. Thus \tf£\k{k)\ 
and \g\ are < 1/i and dephasing is ineffective. In the 
bosonic momentum distribution {hk) was found to relax 
to the one of the corresponding steady state ensemble Qd- 
However, as the dephasing is ineffective, relaxation ^ of 
subsystem density matrices does not occur. This is also 
visible in the observables: As derived in correlators 
do not relax to the value predicted by see also non- 
decaying oscillations of {fix) in Ref. 41 of \S\ and Fig. 10 
of That a particular observable, here {hk), may relax 
anyway is a different issue. In [?! it was shown for a 
slightly modified setup how relaxation of {hk) occurs. 

Discussion. ~ The dephasing theorem ([6|) confirms the 
conjectured (P), clarifies its interpretation, and devises 
conditions for its applicability. Dephasing properties are 
determined in particular by points where the gradient 
g = dkf of the phase function, (fT2|) . vanishes or the 
amplitude /, (fTO|) . diverges. Also note that the notion 
of integrals of motion standing in involution as used for 
classical systems does not carry over to quantum mechan- 
ics, [ic|. Hence it was per se not clear what operators Ik 
were to be chosen in the maximum entropy ensemble ([T|) . 
The theorem settles this question. Further, as ak, Q, 
becomes finite for finite k regions, the subsystem (entan- 
glement) entropies will finally be dominated by the ex- 
tensive contribution Vol log[(l -I- Ce~"'')^'»] (cmp. 
to [m for Id). Hence, the required computational re- 
sources to simulate such systems on classical computers 
scale exponentially in the system size, preventing access 
to arbitrarily long times. On the other hand, this shows 
that quenches are a simple tool for the controlled gener- 
ation of strong (extensive) entanglement. 



Bethe Ansatz integrable systems.- In Bethe Ansatz 
solvable models [12|, the transfer matrix f (A) is conserved 
for any value of the spectral parameter A; [t(A), t(A')] ~ 
and [f(A),-ff] Va,a'- Initial states can be expanded 
in a f (A)-eigenbasis | A) and, via time evolution, nondiag- 
onal contributions will attain quickly oscillating phases 
Qt = EA,A'e(''""''^'^*/'''l^)(^leo|A')(A'|. It will be 
shown elsewhere that, as in the free case, the nondiago- 
nal contributions to the density matrix TrQj_ Qt of a finite 
subsystem £7 will under appropriate preconditions decay. 
Then, the steady state in the thermodynamic limit will, 
in generalization of ([6l), be given by 



limt^oo Trfji Qt = Trj^j_ Qd , Qd 



where p denotes the density of quasiparticles. 

Nonintegrable systems.- Whether or how thermaliza- 
tion occurs in nonintegrable systems is in general un- 
clear. Intuitively, information about the initial state gets 
smeared out by scattering events which are, contrary to 
the integrable case fis'l , able to change the quantum num- 
bers of the involved particles and not factorizable. Our 
results are expected to carry over to nonintegrable cases, 
if system and initial state allow for a description by an 
integrable theory of quasiparticles (e.g. Fermi gases and 
Luttinger liquids) and quasiparticle lifetimes exceed time 
scales necessary to observe dephasing. In such cases, first 
relaxation to the steady state of the integrable theory will 
occur, followed by decay to the thermal ensemble. Nu- 
merical results in [l^ may be interpreted in this vein. 

T. B. thanks the DFG and the Studienstiftung des 
deutschen Volkes for support. 



[1 
[2 

[s: 
[4: 

[s: 
[6: 
[7; 

[9 

[lo; 

[11 

[12 

[13 
[14 

[15 

[16 



E. T. Jaynes, Phys. Rev. 106, 620 (1957). 

T. Kinoshita, T. Wenger, and D. S. Weiss, Nature 440, 

900 (2006). 

M. Rigol et al., Phys. Rev. Lett. 98, 050405 (2007). 

M. Rigol, A. Muramatsu, and M. Olshanii, Phys. Rev. A 

74, 053616 (2006). 

M. A. Cazalilla, Phys. Rev. Lett. 97, 156403 (2006). 

M. Eckstein and M. Kollar, arXiv:0707.2789 

D. M. Gangardt and M. Pustilnik, arXiv:0709.2374 

M. Cramer et al., Phys. Rev. Lett. 100, 030602 (2008). 

A. Varchenko, Funct. Anal. Appl. 10, 175 (1976). 

S. Weigert, Physica D 56, 107 (1992). 

P. Calabrese and J. Cardy, J. Stat. Mech. P04010 (2005). 

N.-C. H. Zachary, Quantum Many-Body Systems in One 

Dimension (World Scientific, Singapore, 1996). 

G. Mussardo, Phys. Rep. 218, 215 (1992). 

C. KoUath, A. M. Lauchli, and E. Altman, Phys. Rev. 

Lett. 98, 180601 (2007). 

One needs to exclude cases where single modes give, in 
the thermodynamic limit, finite contributions to the sum. 
f £ means that J d'^fc] [/(fc)]ij | exists Vij. 



